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Abstract 

In this paper we apply results on the asymptotic zero distribution of the 
Laguerre polynomials to discuss generalizations of the standard large n limit 
in the non-hermitian Penner matrix model. In these generalizations gnU —)■ t, 
but the product QnU is not necessarily fixed to the value of the ’t Hooft 
coupling t. If t > 1 and the limit I = lim„_^oo | sin(7r/5'n)F"' exists, then the 
large n limit is well-defined but depends both on t and on 1. This result 
implies that for t > 1 the standard large n limit with gnti = t fixed is not 
well-defined. The parameter / determines a fine structure of the asymptotic 
eigenvalue support; for I A 0 the support consists of an interval on the 
real axis with charge fraction Q = 1 — 1/t and an /-dependent oval around 
the origin with charge fraction 1/t. For I = 1 these two components meet, 
and for I = 0 the oval collapses to the origin. We also calculate the total 
electrostatic energy S, which turns out to be independent of /, and the free 
energy 4F = S — Qlnl, which does depend of the fine structure parameter 1. 
The existence of large n asymptotic expansions of IF beyond the planar limit 
as well as the double-scaling limit are also discussed. 
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1. Introduction 


The large n limit with hxed’t Hooft coupling ndvr^fyivi = A of non-abelian 
gauge theories has been the subject of intensive research for more than four 
decades, and in particular has fostered the study of large n matrix held 
theories. Random matrix models are a simplihed version of these theories 
that offer a convenient setting to explore large n limits, because the partition 
function of a n x n random matrix model can be written as an integral over 
the matrix eigenvalues x,: 

W{zi) 

/ i=l 

The standard large n limit is dehned here as n —)■ cx) with hxed’t Hooft 
coupling t = ng. This limit exhibits a number of interesting properties in 
the case of Hermitian matrix models with polynomial potentials W(z), such 
as the existence of an asymptotic eigenvalue density supported on a hnite 
number of real intervals (cuts), and the existence of a topological expansion 
of the free energy in the one-cut case PEIEIII]. 

In this paper we discuss this large n limit and generalizations thereof for 
the Penner matrix model with potential 



f (-1 

Jrx-xr V ^ i=i 


W{z) = z + logz. (2) 

The logarithmic function in ([^ is dehned by log 2 ; = In |x| -|- iargx with 0 < 
arg X < 27r and where In denotes the real logarithmic function on the positive 
real axis. We consider the Penner model on any path P in ([^ homotopic in C\ 
[0, -|-cxd) to the steepest descent path \mW{z) = vr through the critical point 
z = —1 of W{z) illustrated in Fig. Obviously, the corresponding partition 
function (Q cannot be interpreted as an integral over the (real) eigenvalues 
of a hermitian matrix model, but to a model in which the integration is 
performed over a set of n x n complex matrices with n eigenvalues constrained 
to lie on the path P. In the literature 0E] these models are frequently called 
holomorphic matrix models. However, to emphasize the relation of ([^ to the 
theory of non-hermitian Laguerre orthogonal polynomials, hereafter we call 
our model the non-hermitian Penner model. 

The hermitian Penner model was introduced in |7j because the large n 
expansion of its free energy provided generating functionals for the Euler 
characteristics Xk,s of the moduli spaces of Riemann surfaces of genus k with 
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Figure 1: Integration contour of the Penner model. 

s > 1 punctures. It was found later j8] that the analytic continuation of the 
free energy in the hermitian case to the non-hermitian case and a suitable 
double scaling limit yielded a generating functional for the Euler character¬ 
istics of unpunctured Riemann surfaces too. As a consequence, the Penner 
model is closely related to the c = 1 noncritical string with one direction of 
the spacetime compactified on the circle with self-dual radius |8l[9l[ini[IIl[l2]. 
More recently, the model has been used to analyze nonperturbative effects 
in c = 1 string theory and matrix models [131 HU- 

In the early 90’s Tan [TillITT] , and Ambjprn, Kristjansen and Makeenko ra 
applied the saddle point method to give what was thought to be a complete 
description of the large n limit of the Penner model for both positive and 
negative fixed values of the’t Hooft coupling. 

The present paper is motivated by the behavior of the asymptotic zero 
distribution of the Laguerre polynomials with negative parameter discovered 
by Kuijlaars and McLaughlin na. Note that hxing the ’t Hooft coupling t 
in effect defines a sequence 

9n = - ( 3 ) 

n 
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which trivially satisfies 


lim uQn = t. (4) 

We will show that the behavior discovered by Kuijlaars and McLaughlin na 
also arises in the large n limit of the non-hermitian Penner model if the 
sequence gn is not restricted to be ([^ , but only required to have the limit Q 

with t > 1 and to have a finite value I for the following limit; 

/ = lim |sin( 7 r/ 5 („)|^/”. (5) 

n^oo 

Under these conditions there indeed exists a large n limit, but this limit 
depends not only on t but also on the parameter /, introduced by Kuijlaars 
and McLaughlin in na. As a consequence, we will also show that the stan¬ 
dard large n limit defined by the sequence ^ with t > 1 does not exist, i.e, 
in this limit there is not an asymptotic eigenvalue density with a well-defined 
support, nor a planar limit of the free energy. 

We begin our analysis by using the saddle point method p Ea [m E] to 
establish the connection between the asymptotic eigenvalue density of the 
Penner model and the asymptotic zero density of the Laguerre polynomials. 
More concretely, we prove that given a sequence of coupling constants Qn, 
the components of the (multiple) saddle points of the integrand of ([^ are 
the zeros of the scaled Laguerre polynomials Therefore 

we can apply the results of na, from which the following picture emerges: 
if f > 1 and the limit I exists, then there exists an asymptotic eigenvalue 
density p{z), and its support 7 ; consists of two pieces: an interval [a,b] on 
the positive real axis and an /-dependent closed loop Ci around the origin. 
These pieces are disjoint for 1^1. If / = 0 the loop reduces to the origin, 
where the asymptotic eigenvalue density has a Dirac 6 contribution. 

In the electrostatic interpretation the eigenvalue density p{z) represents 
a unit-normalized positive charge density on ji in the presence of an external 
electrostatic potential V{z)/t given by 

V{z) = ReW{z), (6) 

and such that the total potential 

U{z) = W{z)-2 f \n\z-z'\p{z')\dz'\ (7) 

is constant on each piece of 7 ^. We calculate explicitly these constants, which 
turn out to be different for / 7 ^ 1. Therefore the density p{z) for / 7 ^ 1 is 
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not an equilibrium density, but a critical density in the sense of Martmez- 
Finkelshtein and Rakhmanov |18) . We also show that the total electrostatic 
energy 

S = - f V{z)p{z)\(\z\ — f |dz| f \(\.z'\\\i\z — z'\p{z)p{z') (8) 


is given by 


Note that 8 is independent of 1. 

We then apply the method of orthogonal polynomials PCD to derive an 
explicit expression of the partition function Q of the non-hermitian Penner 
model in terms of the Barnes G function, and then calculate the planar limit 
of the free energy as 


T = — lim 

n^oo 


1^^ I ^n(^9n) I 




The main result here is that if t > 1 and the limit I exists, then 


( 10 ) 




which involves an /-dependent term. This proves that the I dependence of 
the large n limit also manifests itself in the planar limit of the free energy. 
We illustrate this phenomenon with numerical calculations. 

Finally, by comparing ([^ and ( |II| ) it follows that the difference between 
the free energy and the total energy can be written as 


= £ - (In !)(3[„ i,| 


( 12 ) 


where 


Q[a,b\ 



p{z)\dz 



( 13 ) 


is the amount of charge on the interval [a, b]. This is a somewhat unexpected 
feature of the non-hermitian Penner model; for hermitian matrix models 
with real analytic potentials W{x) that grow faster than In |a:| at inhnity it 
is well known [3j that in the large n limit 


7 = 8. 


(14) 
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We also study in detail the expansion of the free energy beyond the pla¬ 
nar limit for the simple sequence Qn = l/{n/t + a), derive a suitable double 
scaling expansion at the critical value t = 1, and show that this latter ex¬ 
pansion provides generating functions for the Euler characteristics of both 
unpunctured and punctured Riemann surfaces. 

In the last section of the paper we point out some extensions of its main 
resnlts. The technical details of the explicit computation of the total poten¬ 
tial on the support of the eigenvalue density are relegated to an appendix. 


2. Critical eigenvalue densities in the large n limit 

In this section we use the saddle-point method to derive the Schwinger- 
Dyson equation, which is the main tool to determine the asymptotic eigen¬ 
value density p{z), and then discuss the concepts of equilibrium and critical 
densities. The arguments are valid not only for the Tenner model ([^, but 
for any matrix model Q such that W'{z) is a rational function of 2 ; whose 
only singularities are at most a set of simple poles A = {oi,... , 0 ^}- We 
assume that T is a path in the analyticity domain of W{z) such that the 
integral Q is convergent. The Tenner model ([^ corresponds to a case in 
which A = {0}. 


2.1. The saddle point method 

We can write the partition function Q 


as 


^nidn) , 
nl 


„ n 

/ exp ( - .. . , 2 :„) j TTdZi, 

Jrx-.-xr ^ 


where 


•S” = ^ E - 2^ E E 

* * jy* 


(16) 


(16) 


The saddle points z = (zi,..., Zn) of Sn are the solutions of the equations 

(17) 


— W\zi) + '^ —^—=0, i = l,...,n. 

Qn “ Zj - Zi 


Since these equations are symmetric under permutations of the coordinates 
Zi, generically each solution z = {zi,..., Zn) of 0 gives rise to a set of 
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n\ solutions ct(z) = {za(i), ..., ^cr(n)), where a denote the permutations of n 
objects. 

The saddle point method assumes that there exists a sequence of saddle 
points of Sn and a unit-normalized positive density p{z) 

(the eigenvalue density) with support 7 such that 


p{z)\dz\, 


n 


( 18 ) 


2=1 


in the large n limit (|^. Each saddle point of the sequence determines a 
monic polynomial which we denote by 


Sn ( z ) = P [(2 - zf ’). 


( 19 ) 


and, as a consequence of 0 , the resolvent function 

= - y" —^ 

n Sn{z) z- 


satisfies the Ricatti equation 

icz'(z) +a.„(z)" - —W'{z)^,.(z) ^-4-11 

n npn n^Pn ^ 


W'{z) - IT'(4”^) 


z — z. 


(n) 


( 20 ) 


( 21 ) 


This equation coincides with the Ricatti equation (2.6) derived in |17( after 
performing the notational replacements N ^ n, g ^ ng^ and u —)■ — 


The large n limit oj{z) of the resolvent 0 Jn{z) defined in (20) is 

p{z')\(lz'\ 


Uj{z) = 


z — z' 


( 22 ) 


Therefore, as a consequence of (21), the function u{z) must satisfy the 
Schwinger-Dyson equation 


uj{zy — ^W' {z)oj{z) = —- 


W'{z) - W'{z') 


t 


t 


z — z' 


p{z')\dz'\. 


(23) 


For later reference, we also show the reduction to a second order linear equa- 

Sn{z). 


tion for Sn{z) of the Ricatti equation (21); 

W'{z) - IT'(4”^) 


S^z)-j-»"(z)Si(z) = -j-(5; 

Pn Pn ' 


2 : — 7- 


(n) 


( 24 ) 
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2.2. Critical densities 

In this section we discuss briefly the minimum technical background that 
permits us to differentiate between the equilibrium density (i.e, the absolute 
minimum) and critical densities in the sense of Martmez-Finkelshtein and 
Rakhmanov [H]. In this context, it is convenient to introduce the function 


y(z) = -W'{z) - 2cj(z), 


and rewrite the Schwinger-Dyson equation (23) as 

p(z)^ = R(z), 


(25) 


(26) 


where 


= (^W'(z)) - - 


4 fW'{z)-W'{z') 


z — z' 


p{z')\dz'\. 


(27) 


Incidentally, the density p{z) can be recovered from y{z) using the Sokhotskii- 
Plemelj formulas. Moreover, under our assumptions on W'{z) the function 
R{z) defined in (27) is a rational function of ^ with the same poles as W'{z). 
Hence, from (25) it follows that the function y{z) is analytic outside 7 U M, 
and (26) implies 


!/(2+) = -!/(2-), 2 e 7- 


(28) 


Equation (26) is crucial to derive the main features of the eigenvalue 
density. Thus, we note that the function u}{z) coincides with the Cauchy 
transform 

dp{z') 


Jc ^ ^ 

of the measure on C, with support 7 , given by 

dp{z) = p{z) |d^|. 


Hence (26) can be written as 

T 


{pv'(z) - 2C‘‘(z)) = R(z). 


(29) 


(30) 


(31) 


This relation shows that dp{z) is a continuous critical measure on C in the 
sense of Martinez-Finkelshtein and Rakhmanov [18]. As a consequence (see 
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Lemma 5.2 of [TS] and Proposition 3.8 of M). the snpport 7 of p{z) is a 
nnion of a finite nnmber of analytic arcs 


7 = 7 i U72 U • • • U 7 s, 

which are maximal trajectories of the qnadratic differential 

-y{zf{dzf = -R{z){dz)\ 
i.e., maximal cnrves [ 20 ] z = z{t) (t G {a,/3)) snch that 

>0, for all ( 6 (a,/3). 
Moreover IISl II HI . the total potential 


U(z) = p'(z) - 2 


In \z — e'|p(o')|do'|, 


(32) 


(33) 


(34) 


(36) 


(36) 


(37) 


is locally constant on 7 

U{z)=Ui, z G ji, i = 

with possibly different constants Ui, and that it satisfies 

dU{z) dU{z) 
dn+ ~ dn. ’ 

where n± denote the two normal vectors to 7 at z pointing in the opposite 
directions. Eqnation (37) is the so-called S'-property of Stahl PH 1221 [23] 
and of Gonchar and Rakhmanov PIES], whose electrostatic interpretation 
is that the electric helds at either side of 7 are opposite or, eqnivalently, that 
the forces acting on each element of charge at z from the two sides of 7 are 
eqnal. 

The total electrostatic energy ([^ can be written in terms of the external 
potential and the constants Ui as 


£ = — I V{z)p{z)\dz\ + - / U{z)p{z)\dz\ 


1 

1 


V{z)p{z)\dz\ +-^Ui / p{z)\dz\. 


(38) 


2=1 


'li 
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We emphasize that unless all the constants Ui are equal, the eigenvalue 
density is not an equilibrium density minimizing the total energy ([^. In 
fact, continuous critical measures are characterized by a different stationary 
condition for the total energy, 

DhS[p] = lim = 0 , (39) 

s 

under variations p{z) —)■ p^{z), with d/i^( 2 ;*) = dp{z), induced by local set 
transformations z^ = z sh{z) induced by smooth functions h{z). 

Summing up, in general the eigenvalue densities arising from the saddle 
point method determine continuous critical measures on C, and we will refer 
to them as critical densities. 


3. Laguerre polynomials and the asymptotic eigenvalue density of 
the Penner model 


In this section we specialize the general equations obtained in the previ¬ 
ous section to the Penner model ([^ and show its relation to the theory of 
Laguerre polynomials [ 8 ] 




k=0 


n + a \ {—z)^ 


n — k 


k\ 


Indeed, the saddle point equations ([T^ for the Penner model 

2 

= 0, ? = l,...,n. 


are 


1 f 1 
-(l + 

9n V 


(n) 


Z- 


E 


M) (n) 




and the corresponding Ricatti equation (21) is 


—Uj'{z)+Un{zy - (l + -^Uniz) — — ^ ... 

n ngn\ zJ n^9nZ^zl’ 


E 


From (41) it follows that 


^ = -n, 


i=l Zi 


(n) 


and we get the following second order linear equation (24) for Sn{z): 

1(1 + 1 

9n ^ Z^ 


- l(i + 

CJn ' Z ' j/?2 ^ 


(40) 

(41) 

(42) 

(43) 

(44) 
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By comparing (44) with the Laguerre differential equation 


zu''{z) + (a + 1 — z)v!{z) + nu{z) = 0 , u{z) = L^^\z), 


( 46 ) 


we find that the monic polynomials Sn{z) are proportional to the rescaled 
Laguerre polynomials 


V( 7 „/ 


CX-n 1 

9n 


(46) 


Therefore, the saddle points of the Penner model with 

coupling constants gn are given by 


An) 


9n 

where are the zeros of L^n\z) 


{an,n) 


i = 1,... ,n, 


(47) 


3.1. Zero asymptotics of scaled Laguerre polynomials 

Riemann-Hilbert and steepest-descent methods Ha EHl EH EH] have per¬ 
mitted the complete characterization of the asymptotic zero distribution 
Pl{.z) of the scaled Laguerre polynomials Llf‘"\nz) for all real values of the 
parameter 

rv 

A = lim —. (48) 

n—^oo Tl 

We will concentrate on the case — 1 < ^4 < 0. It turns out that the zeros 
cluster along certain curves 7 , but the form of 7 depends not only on the 
value of A but also on the parameter 

1= lim I sin(a„ 7 r)|^/” = lim [dist(a„, (49) 

n^oo n^oo 

which represents the proximity degree of the sequence an to the integers. 
Thus, if the limit I exists and we denote 

a± = A + 2± 2^41 + 1, (50) 


then (see Theorem 1.2 and Eq. (3.14) of [27]) it follows that (x[a_,a+](a^) 
equals 1 if a; G [a_, a+] and zero otherwise): 

1. For I = 0 

7 = {0} U [a_,a+], (51) 

and 

PL{z)\dz\ = -A6{z) + «-)(«+ (52) 
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2. For I > 0 


7 = F; U [a_,a+], (53) 

where F^ C C \ ({0} U [a_, +C)o)) is a simple closed curve encircling 0 
clockwise, which is determined by the implicit equation 


Re 




dz' = - In/. 


(54) 


Furthermore 


Pl{z)\<1z\ = dpi{z) + X[a-,a+]{.x)dx, (55) 


2t^x 


where 


, / X \/{.z-a_){z-a+) 

d/x«(^) = -7777- 


27ri7 


(56) 


and the branch of the squared root is such that \/[z — a-){z — a+) ~ 7 
as 7 —)■ cx). 

The case / = 1 is generic (see remark 1.3 in jl5| ) and it is the only value 
of / for which the loop F; and the interval [a_, a+] intersect (at the point a_). 

We remark that this fine structure of the asymptotic zero distribution of 
the Laguerre polynomials does not arise for other values of the limit (48). 
Thus, if —cx) < A < —1 the support of the zero distribution is a well-defined 
simple analytic curve symmetric with respect to the real axis 

3.2. The asymptotic eigenvalue density of the Penner model 


In view of (46) and (47), by setting 
1 


i = -^. 71 = {*2, 2 £ 7 }, = 


(57) 


we can apply immediately the results on the asymptotic zero distribution of 
Laguerre polynomials into properties for the asymptotic eigenvalue density 
of the Penner model. 

Thus, it follows from (53)-(56) that for f > 1 and I > 0 the support 'ji of 
the eigenvalue density is given by 


= G U [a, b], 


(58) 
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where 


a = 2t-l- 2Jt{t-l), b= (59) 

a 

and Cl G C\ ({0} U [a,+oo)) is the simple closed curve determined by the 
implicit equation 


Re 

Furthermore, the eigenvalue density is 

p{z)\dz\ = dpi{z) + (61) 

where 

with 

r{z) = \/(z — a)(z — b) (63) 

such that r(z) ^ z as z ^ oo. 

We will now show an important difference of these large n limits with 
respect to the hermitian case: for I ^ 1 the eigenvalue density is not an 
equilibrium density but a critical density. Indeed, let us write the total 
potential in the form 

U{z) = -V{z)-2Reg{z), (64) 

where g{z) is the function |15) 




(60) 


9 (z) = / log(2 - 2')p(z')|d2'|. 


(65) 


'll 


For z' G [a, b] we take the cut of log(^ — z') along [z', cx)), while for z' G Ci 
we take the cut along the part of Ci that starts at z', continues along Ci in 
the clockwise sense until it reaches the real axis and then runs to +cx). Thus 
the function g{z) is analytic in C \ {Ci U [a, cx))). 

Let Doo denote the unbounded connected component of C \ {Ci U [a, cx))). 


Then from (61)-(63) and using Cauchy’s integration it follows that 




z G T)qq , 


( 66 ) 
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and as a consequence we have that for z G D^o 


Re qiz) = — fx + In R| — a — Ina^-Re 

' 2t\ ' ' ) 2t 


+ Re5f(a). (67) 


Hence from (60) and (64) we obtain 


U{z) = 


-(a + Ina) — 2Re5f(a) — In/, for z G 7 ;, 


( 68 ) 


-(a + Ina) — 2 Re 5 f(a), forzG[a,6]. 
^ t 


Therefore 

= (69) 

Since for I ^ 1 the constant values of the potential on each component of 
the support are different, the eigenvalue density p{z) is not an equilibrium 
density but a continuous critical density. We defer to appendix A the explicit 
evaluation of the constant value on the real interval, which turns out to be 

U\[a,b\ = (2 - - In/ - (1 - ln(/ - 1). (70) 

3.3. Existence of the eigenvalue density in the large n limit 

It is worth noticing several striking consequences that the above results 
imply for the large n limit (|^ of the Renner model. We have seen that if 
/ > 1 and the limit I exists, then the asymptotic eigenvalue density of the 
Renner model exists but depends not only on the value of the parameter t 
but also on the value of 1. We will prove now that in the large n limit with 
hxed’t Hooft coupling / > 1 there are subsequences of the sequence gn = t/n 
that lead to different values of I and therefore determine different asymptotic 
eigenvalue densities. Consequently, the large n limit of the eigenvalue density 
corresponding to the sequence gn = t/n does not exist. 

We hrst observe that in terms of the sequence of coupling constants the 
parameter / can be written as 


/ 


lim 

n^oo 



1/n 


lim 

72^00 




( 71 ) 


where {x} stands for the fractional part of the real number x. 
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Let us suppose that the fixed ’t Hooft coupling t takes a rational value 
(in lowest terms) 

t = -, {p> q). (72) 

Q 

Then nq = jVi(modp) with jVt = 0, — 1, and in the sequence of 

fractional parts {1/gn} = each number 0,1/p,..., (p — l)/p, appears 

infinitely often. Therefore there are subsequences of {l/p„} with I = 0 
and subsequences with 1 = 1. Likewise, if t is irrational then the sequence 
{1/gn} = {n/t} is a dense subset of the interval [0,1]. Hence {l/p„} has 
subsequences with I = 0 and subsequences with / = 1. 

There are simple examples of large n limits (|^ such that the limit I exists. 
For instance, take t rational of the form (72) and define 


gn — 


n/t + a ’ 


a > 0. 


(73) 


Decomposing again nq = jn (modp) with = 0,1,... ,p — 1, it follows that 
the sequence {1/gn} takes the values {a} , [a + 1/p} ,... {a + (p — l)/p}. 
Hence for real values of a > 0 such that 


k 

a ^ K = 1, 2,' 

p 


(74) 


the sequence yields the value / = 1. 


Unfortunately, the example (73) does not work for t irrational. However, 


we may provide a more general example. Take the large n limit with 

1 e-™ 


g-n 


[n/t] + Cn 


Crt — 


n = 1,2,... 


(75) 


where t > 1 and [x] denotes the integer part of x. We assume 0 < r < +cx) 
and e~°° = 0. Then it is obvious that p„n —)■ t and that 


/ = lim = 

n^oo 


(76) 


Therefore, all the possibilities (|51[)-(|56[) arise. 

We may use this latter example to illustrate numerically the dependence 
of the support on 1. Figure [^corresponds to a sequence (75) with t = a/3 and 
r = 1/7. The red markers denote the zeros of the corresponding scaled La- 
guerre polynomial with n = 60, which already accumulate distinctly on the 
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Figure 2: Zeros of the n = 60 scal ed Laguerre polynomial and asymptotic density support 
for a sequence of the form (751 with t = -s/S and r = 1/7. 


two components Ci and [a, b] of the support, which have been calculated nu¬ 
merically from (59) and ( [60| . Figure [^corresponds to a similar sequence (75) 
with t = y/S but r = 1/3. Note how the endpoints (59) a ~ 0.21 and b ~ 4.72 


(marked in blue) of the component of the support on the real axis are the 
same in both instances, but the oval Ci is notoriously different (both figures 
are drawn with the same scale). Note also how in both cases 25 out of the 
60 zeros are on the real interval, i.e., a fraction 25/60 ~ 0.417 already quite 
close to the limiting value 1 — 1/t = 1 — 1/a/3 ~ 0.423. 

Finally, we notice that the results of ra on the asymptotic zero distribu¬ 
tion of Laguerre polynomials for values —oo < A < — 1 of the limit (48) can 


also be applied to the non-hermitian Fenner model. These results imply the 
existence of an asymptotic eigenvalue density with a well-defined support for 
any sequence gn of coupling constants with limit (j^ such that 0 < f < 1. 

4. The large n limit of the free energy 

The method of orthogonal polynomials PE] 


J Pn{z)z^e = 0, /c = 0,..., n — 1, 

leads to the following expression of the partition function 


n—1 




,n—k 
k 1 


(77) 


(78) 


k=l 
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0.10 


t=yf3, ^ 1/3 


2 3 4 

Re z 


Figure 3: Zeros of the n = 60 scaled La guerre polynomial and asymptotic density sup¬ 
ported for a sequence of the form (751 with t = VS and r = 1/3. 


where 


ho = J (79) 

and Tk are the recurrence coefficients in the three-term recursion relation 

zPk{z) = Pk+i{z) + SkPk{z) + rkPk-i{z). (80) 

For a < — 1 the Laguerre polynomials satisfy an orthogonality relation of 
the form daEZ] 


Ll^\z)z^z°‘e ^dz = 0, fc = 0,1,..., n — 1, 


(81) 


where F is a path of the form shown in Fig. and 2 ;“ = exp(Q! log z) with the 
branch of log defined in the Introduction. It is immediate that the rescaled 
Laguerre polynomials 

4-'/»)(^) (82) 


are proportional to the orthogonal polynomials Pn{z) (77) corresponding to 
the Fenner model We apply ([^-([^ to calculate the partition 

function of the Fenner model. 

Thus, we have that 


hn = / 2 ; se adz. 


(83) 
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This integral is an analytic function of 5 ^ in C \ {0}. For Keg > 1 we may 
evaluate ho by approaching the contour F to the positive real line and then 
extend the result analytically. Then, for positive values of g we get 

/io = ( 1 (84) 

Moreover, from the three-term recursion relation of the Laguerre polyno¬ 
mials 

= -{n + + {2n +a + 1 )LI^\z) - (n + a)L^^l^{z), (85) 

the recurrence coefficients are given by 


rk = kg{kg - 1 ). 

At this point it is convenient to use the Barnes G function 
by the canonical product 


G{z + 1) = (27r)^/V5("+^'(i+^)) Yl{l + ^ye 


k=l 


^—z+z^l2k 


where 7 is the Euler-Mascheroni constant. This function satisfies 
G{z + l) = T{z)G{z), G(l) = l, 

and therefore 


n— 1 


ll(k + 


a 


^n-k -F a -f-1) 


k=l 


G(a + l)F(a + l)' 


( 86 ) 

defined 

(87) 

( 88 ) 

(89) 


Substituting (84) and ( 86 ) into ( |78[ ) and using (89) we find that the partition 
function can be written as 


Zn(9)=9”‘”-'’(l 


e 3 


G{n + l)G(n -|- 1 — 7j 


G(l-7 

9 


(90) 


The second factor in this expression is an essential feature of the class of 
non-hermitian Fenner models. In order to analyze the large n limit (|^ we 
write it as 


l-e-‘T) = e“T"(2i)”sin’^ ( - 


TT 
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(91) 







As to the third factor in (90), when we evaluate Zn{gn) for a sequence 
Qn we need to evaluate the Barnes G functions for large positive and large 
(in absolute value) negative arguments. More concretely, if f > 1 the terms 
G(n +1) and G{n + 1 — 1/ Qn) have to be evaluated for large positive values of 
their respective arguments, while the term G(1 — 1/gn) in the denominator 
has to be evaluated for large (in absolute value) negative values. However, 
if 0 < f < 1 only the term G{n + 1) has to be evaluated for large positive 
values, while both G{n + 1 — 1/gn) and G(1 — 1/gn) have to be evaluated for 
large (in absolute value) negative values. 

The asymptotic expansion of G{x) for large x > 0 is well-known and often 
written keeping an unexpanded gamma function (see equation 5.17.5 in |3U)). 
but we hnd more convenient the fully expanded version 


1 3 1 1 

logG(x + l) ~ -x^logx--x^ + -xlog(27r)-—logx + C'(-l) 


E 

k=2 


B 


2k 


2k{2k-2) 


X —>■ oo. 


(92) 


where i? 2 fc are the Bernoulli numbers and is the derivative of the Riemann 
zeta function. Incidentally, the numbers Xk,o = B 2 k/{‘^k{2k — 2)) which 


appear in (92) are precisely the virtual Euler characteristics for the moduli 

To relate the values of the 


space of unpunctured Riemann surfaces 
Barnes G function on the negative axis to the values on the positive axis we 
use the following expression m, 


G{-x) = (-l)L^/2J-i 


G{x + 2) 


Sm TTX 


TT 


x+1 


g^Cl2(27r(3;-[a:J)) 


X > 0, (93) 


where Cl 2 (x) is the Clausen function (called Clausen’s Integral and denoted 
by / in |33]). Therefore, for sequences gn such that the limit ([^ exists, taking 
into account that the Clausen function is bounded, and that 


In gn 
gn 


= - In 
t 


-) +o(n), 
n/ 


n —>■ oo. 


(94) 


it follows from (90)-(93) that the free energy (10) is 

in-1 




, , 1, 3 /f- 1 

In /-In f H— (- 

2 2\ t 


)4(4^) H-V (-) 


where H (x) is the Heaviside step function. 
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Figure 4: Numerical values of for a sequence with t = -s/S and r = 1/3 in (751 (upper 
one, blue markers) and a similar sequence with t = -s/S and r = 1/7 (lower one, red 
markers). The horizontal lines are the respective limits given by (95). 


As a numerical illustration of this result, in Fig. |^we show the values of 

1^^ I i.9n} I 


T = 

n. 




(96) 


for two sequences Qn- the upper one (blue markers) corresponds to t = a/3 
and r = 1/3 in (75), while the lower one (red markers) corresponds to f = Vs 
and r = 1/7. In this example f > 1 and the horizontal lines are the respective 
(different) limiting values given by (95). 

Similarly, in Fig. [^the upper, blue markers correspond to t = 1/a/3 and 
r = 1/3 in (75), while the lower, red markers correspond to t = 1/a/3 and 
r = 1/7. In this case t < 1 and the horizontal line is the common limit given 
by 


4.I. Free energy versus total energy 

For f > 1 the expression (95) for the planar limit of the free energy reads 

1 /f - 1' 


tF = [ - - --hit 


We are going to prove that this identity can be written as 


(97) 




(98) 
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Figure 5: Numerical values of Tn for a sequence with t = l/-\/3 and r = 1/3 in (75) (upper 
one, blue markers) and a similar sequence with t = I/'s/s and r = Xjl (lower one, red 
markers). The horizontal line is the common limit given by 


where Q[a,b] is the charge contained in the interval [a,b], 


Q[a,b] = / p{z)\dzl 


and £ is the total energy. 


Using Cauchy’s theorem and (62) we obtain 

r(0) 


p{z)\dz\ = 


'Cl 


Hence 


p{z)\dz\ = 1 - 


In this way the first term in ( |97| ) can be expressed as 

(In /) Q[a,b] ■ 


From (38) and (69) we have that the total energy is 

1 
2 


^ ^ / V{z)p{z)\dz\ + 




(99) 


( 100 ) 

( 101 ) 

( 102 ) 

(103) 


The integral term in (103) can be also explicitly calculated (see appendix A) 

. (104) 


'll 


V{z)p{z)\dz\ = In/ + t — 2 + ^1 — -j ln(t — 1), 
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and the final expression for the total energy takes the form 

1 . 


^ . 3/t-l 




( 105 ) 


Therefore, the identity (98) follows. It mnst be noticed that, nnlike the free 


energy, the total energy is independent on the valne of 1. 

4-2. Beyond the planar limit. The double scaling limit 

It has been rigoronsly proved for a large class of hermitian models with 
polynomial potentials BElll] that in the one-cnt case the partition fnnction 
admits an asymptotic large n expansion of the form 


— In 






kU 


2-2k 


(106) 


fc =0 


where denotes the partition fnnction of the Ganssian model W[z) = 12. 

It is also well-known the appearance of oscillatory terms in the mnlti-cnt 
case |M|. The qnestion natnrally arises if similar expansions can be generated 


for the non-hermitian Tenner model from (90) and snitable seqnences Pn 
satisfying (|^. 

Let ns consider the one-parameter family of seqnence of conpling con¬ 
stants 

a > 0. (107) 


9n — 


n/t + a 


As we have seen above in Sec. 3.3, if f is a rational nnmber t = p/q > 1 
and a ^ k/p {k = 0,1, 2,...) then the seqnence pn determines a well defined 
large n limit with / = 1. Snbstitnting this seqnence pn in the exact eqnation 


for the partition fnnction (90) and expanding in n we obtain 


— \n\Zr, 


— In n 

12 


E 

fc =0 


n 


2-k 


where the coefficients are given by: 

1. 3/f-l 

Fq = — In f H— (- 

“ 2 2\ t 


Fi = — ln(27r) — a [ 2 — t — 


t - 1 


t 


ln{t - 1), 


ln(f — 1) 


(108) 


(109) 


t - 1 


In 


2sin I TT 1-ha 

p 


Jr, 


( 110 ) 
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F, = -C'(-l) + ^M*-l)-y(«(*+l)+ln(*-l)) 


a In 


2 sin ( TT ( — + a 
P 




F, = 


.l)fc+l^A:^fc-l(fc 1)^+1 


Bk 


k{k — 1) k{k — 2) 


,k-2 



1 ) 


k-2 


-1)^ 


2i 


-2j(2j-2) 


fc - 3 
2j-3 


a 


X 


fc-2j 


a 


+ 


a 


k-2 


k{k-l){k-2) l2{k-2) 

{k > 3), (112) 

and jn = nq {mod p). Note that the expression for Fq coincides with the 


F of (95) for f > 1 and I = 1. Note also that this expansion (108) features 


oscillatory terms in the coefficients Fi and F 2 . The coefficients F^ with k >2 
diverge as f —)■ 1+, so that the limit value f = 1 is critical. 

In order to study the free energy near the critical value t = 1 and to avoid 


the limitation of using rational values of f > 1, we consider the sequence (107) 


for 0 < f < 1 and arbitrary values of a > 0. In this case, substituting the 


sequence Pn in (90) and expanding in n we obtain an asymptotic expansion 


of the form (108) with the same coefficients Fk for k >3 but 

1 
2 


r. 1. 3/f-l' 

“2 2\ t 




Fi = — ln(27r) — a ( 2 — f — 


t 


ln(l — t) 


t 

F 2 = -C'(-l) + — ln(l -t) - — {t{t + 1) + ln(l - t)) 


(113) 

(114) 

(115) 


Note that the oscillatory terms arising in the expansion for f > 1 have dis¬ 
appeared. At f = 1 the coefficient F 2 has a logarithmic singularity while the 
coefficients Fk with k > 3 have a pole of increasing order: 


E B2j 


k — 3 


(^-l)"“n^2j(2j-2) V2J-3 


a 


k-2j 


a 


+ 


a 


k-2 


k{k-l){k-2) I2{k-2)J' 


k>3. 


( 116 ) 
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The sum of the divergent terms with k > 3 (as t —)■ 1+) of the asymptotic 
expansion of the free energy can be formally written as 


In Zr® 




k=3 

oo 


k-2j 


- 1 )) 


2-k 


k=3 


a 


a 


k-2 


k{k-l){k-2) 12{k-2)J' 

(117) 


It is then natural to regularize this expansion with a double scaling limit in 
which f —)■ 1 in the interval 0 < f < 1 and n ^ oo with 


{t — l)n = ^ = fixed. (US) 

We sum separately the terms with even k = 2h and odd k = 2h + l in the first 
series in (117) by changing sums in (fc, j) into sums in (s, j) with s = h — j. 
The second series in (117) can be evaluated as a sum of iterated integrals of 

(119) 


the geometric series. Finally, in terms of 


a 


we find that 


where 


In _ /(;x,r) + 

i=2 


( 120 ) 


4 


and 


with 


s=0 


B 


-r)2ln(l-r)) + ^ln(l-r) 

(121) 

;(-l)‘+‘x',,.r', 

) 

(122) 

' 2 j / 2 j — 3 + s \ 

'-2)V 5 ;■ 

(123) 


Xj,s = (-1)' 


These Xj,s are precisely the Euler characteristics of the moduli space of Rie- 
mann surfaces of genus j > 2 and s > 0 punctures. Therefore the double 
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scaling expansion (120) comprises generating fnnctions of the Enler charac¬ 


teristics not only for the unpnnctured Riemann snrfaces, bnt also for the 
punctnred ones. This featnre is a conseqnence of the dependence of the 
sequence Qn on the parameter a. In particular, if we set r = 0 we find 


- 1 " Fr‘L.„ ~ 


i=2 


2-2j -^2j 


2j(2j-2)’ 


(124) 


which reproduces the term of genus j > 2 of the topological expansion for 
the free energy of the c = 1 string. 

5. Concluding remarks 

In this paper we have discussed generalized large n limits for the Penner 
model, and have considered the value t = 1 only through a limit t —>■ 1 of a 
particular example. However, the case t = 1 can be independently studied 
in terms of the asymptotic zero distribution of Laguerre polynomials in the 

T. This zero distribution was determined in [28] and 
(see Theorem 1 of [28]). Indeed, the support 

but now through the 


limit (48) with A = 


also exhibits a fine structure 
of the zero density p{z) depends on the sequence a. 
parameter 

m = lim [dist(a„,Zn)]^/”, 


(125) 


where Z„ = { — 1, —2,..., —n}. As a consequence, there is also in this case 
a one-parameter family of large n limits of the Penner model for t = 1 


depending on the value of m. However, since in general the parameters I (49) 


and m (125) are different (take for instance the sequence «„ = — (n -1- 1) -|- 
(2/3)"'), the large n limit of the Penner model for t = 1 requires a separate 
analysis. 

Multi-Penner models of the form 


- 9i). 


(126) 


2=1 


are relevant to characterize the correlation functions of the d = 2 conformal 
Ai Toda field theory Therefore, it is interesting to investigate their large 
n limits [36] • The simplest case is 


W(z) = Pi log(z -1)+ P 2 log(^ 1), 


(127) 
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which is closely related to the theory of Jacobi polynomials 

'n + 13' 




k=0 


n + a 
n — k 


k 




n—k 


(128) 


The asymptotic zero distribution of Jacobi polynomials in the large n 
limit 

n —)■ oo, — —)■ A, — —)■ S, (129) 

n n 

for has been determined in ra and |3H]. It turns out that, as in 

the case of Laguerre polynomials, the zero density and its support depend 
not only on A and B but also on additional parameters related to the degree 
of approximation of the sequences and f3n to the set Z of integers. In 


this way the large n limit of the two-Penner matrix model (127) exhibits a 


fine structure similar to the Penner model studied in the present work. We 
finally mention that the analysis of these features for general multi-Penner 


models (126) would require the characterization of the asymptotic zero dis¬ 


tribution of Stieltjes polynomials |18) . 


Appendix A 

To compute Re 5^(0) as a function of t we write 

I C . , r(z) 


Reg (a) = Re 


27rit 


log(a 


dz 


'Cl 


, 1 f 

H-/ ln(x —a)-^^- dx. 

2TTt Ja ^ ’ X 

The first term reduces to (In a) / 1. As for the second term 
"TT J^ 


X 


^la 


6 -|- 2 In -h (a -|- 6) In — 


Taking into account that b = Ija^ from (68)-(131) we get 


Reg(a) = — 


a -2 In 2(0-1— ) — fa-|- 2 ) In a 

a \ a j \ a 


(130) 


(131) 


J- -I -1-2^ ln(l J- a) J- -I -2^ ln(l — a) 


,(132) 
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so that 


U\[a,b] = -{a+ \n a) -2Re g{a) 


t 

1 


CL J (l-|-2 log 2) -j- ^ci -j- — J In CL 


CL -\ -h 2 

a 


j ln(l + a) 




a H-2 

a 


1 ln(l 




Thus, using the relations 

a H— = 4t — 2, 
a 

2 ln(l + a) — In a = In t + 2 In 2 , 

2 ln(l — a) — In a = ln(t — 1 ) + 2 In 2 , 
it follows at once that 

£'lMl = (2-7)-l"i-(l-7) Mi-1). 

Let us now consider 

f V(z)p(z)\Az\ 


(133) 

. (134) 

(135) 

(136) 

(137) 

(138) 


'7i 


'll 


In | 2 ;|p(^)|d 2 ;|-— / \/{x — a){b — x)da; 
27rt L 


= Reg{0) + 


{b — of 


m 

= Re ( 7 ( 0 ) - 1-1 — 1. 


(139) 

(140) 

(141) 


The computation of Re 77 ( 0 ) requires a more delicate treatment. We observe 
that 


Re( 7 (^) = / \yi\z — z'\p{z')\dz'\ 


(142) 


'll 


is a continuous function on C which is harmonic in C \ 7 ^. We already know 
that (67) is the form of this function in the unbounded connected component 
Hoo of C \ {Cl U [a, cxd)). To determine Re( 7 (^) in the bounded connected 
component ZIq of C \ {Ci U [a, 00 )) we introduce the function 


t( -I- 1 

f{^)— / -d 2 ;' + lna, z&C\[a,oo), 


(143) 
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which is analytic on the complex plane with a cut along [a, oo). We have 
that 


In | 2 ;| = tlnZ + Re/( 2 ;), z E Ci. 


Then from (67) we obtain 


Reg{z) = — — a — In a + Re f{z)^ + Re 5 f(a) + In /, z E ■ji. 


(144) 

(145) 


Since the right-hand side of this equation is a harmonic function in C\ [a, cxd), 
it follows that the same expression for Re 5 '( 2 ;) holds in Dq. As a consequence 
we find 


Re 5 f( 0 ) = — ^ — a — Ina-|-Re/(0)j-|-Re 5 f(a)-|-In/, 


Moreover, we have that 

Re/(0) = lna-|- 


^(a — x)(& — x) — 1 


dx 


= 21n2 — l-|-lna-|- 




X 

1 / 1 
9 + “ 

L2 V a 


ln(l — a) 


ln(l -|- a) 
= —1-|-(t — 1) ln(t — 1) — t Int. 


Substituting (138) and (148) into (146) then (139) implies that 

1 ' 


'71 


V{z)p{z)\dz\ = Inl + t — 2 + ^1 — - j ln(t — 1). 


Hence, from (103), (138) and (151) we finally get 

^ 1, 3/f-l\ l/f-l\2 

£ = —lnt+- - -- - 

2 2\ t J 2\ t 


ln(t — 1). 


(146) 

(147) 


(148) 


(149) 

(150) 


(151) 


(152) 
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